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ABSTRACT 


A  modal  has  been,  developed  for  calculating  x-ray  line  emission 
from  spherical  plasmas. 

Hie  main  features  of  this  method  are: 

1)  Plasma  parameters  are  obtained  from  a  one-dimensional 
Lagrangian  hydrodynamics  and  heat  flow  code. 

2)  Multi-frequency  groups:  the  line  structure  can  be  reproduced 
with  the  desired  accuracy  by  adjusting  the  number  of  frequency  groups. 

3)  Self  consistent,  time  dependent  excited  level  populations  and 
radiation  fluxes:  the  code  starts  with  coronal  populations,  calculates 
the  ensuing  radiation  flux  and  then  recalculates  the  populations  and  so 
on,  iterating  until  convergence  is  reached. 

4)  Geometrical  groups  of  rays  grouped  by  spherical  impact 
parameter. 

3)  Line  broadening  due  to  ionic  thermal  agitation  and  Doppler 
shift  due  to  the  net  plasma  flow  velocity.  Inclusion  of  the  flow 
velocity  shift  vould  be  difficult  without  the  multi-frequency  group 
treatment. 

The  method  has  been  applied  to  an  aluminum  target,  and  the 
results  are  in  good  agreement  with  previous  experimental  work.  The 
total  energy,  summed  over  all  lines,  as  well  as  the  line  intensity 
ratios  (which  are  a  sensitive  measure  of  agreement  with  experiment)  were 
predicted  with  good  accuracy.  The  pictures  that  would  be  seen  by  a  pin¬ 
hole  camera  are  also  calculated  by  the  code. 
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CHAPTER  1 


INTRODUCTION 

The  scope  o £  this  work  is  to  calculate  the  transport  in  and 
emission  of  x-rays  from  a  spherical  plasma  created  by  laser-heating  of  a 
spherical  target. 

The  setting  for  such  experiments  is  as  follows:  a  spherical 
target  is  hit  from  every  direction  by  a  laser  pulse  split  into  a  number 
of  beams  so  as  to  approximate  spherical  symmetry.  Due  to  the  intense 
energy  deposition,  part  of  the  material  is  instantly  vaporized  and 
ionized,  thereby  creating  a  plasma  cloud  of  density  decreasing  with 
radius.  As  is  known,  light  of  frequency  f  cannot  penetrate  a  plasma  of 
electronic  density  equal  to  or  greater  than  a  'critical  density'  at 
which  the  plasma  frequency  equals  the  frequency  of  the  incoming  light. 
The  radius  at  which  the  electronic  density  in  the  spherical  target 
plasma  equals  the  critical  density  for  the  frequency  of  the  laser  light 
defines  the  'critical  surface.'  Energy  is  transported  inward  of  the 
critical  surface  mainly  by  electronic  thermal  conduction.  The  intrinsic 
non-linearity  of  the  heat  transfer  in  plasmas  produces  a  thermal  front 
propagating  inward:  the  thermal  wave.  The  temperature  distribution  has 
therefore  two  opposite  faces:  from  the  thermal  front  out  the  high 
conductivity  of  the  electrons  makes  the  spatial  behavior  of  the 
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temperature  very  smooth,  vheren  inward  of  the  thermal  front  the 
temperature  drops  abruptly  and  the  plasma  is  essentially  cold. 

xhe  hot  part  of  the  plasma  is  a  very  strong  x-ray  emitter,  and 
the  present  work  has  as  goal  to  predict  what  an  observer  looking  at  the 
plasaa  from  outside  would  see. 

This  task  is  two-fold:  modeling  the  hydrodynamics  of  the  process 
and  modeling  the  x-ray  transfer.  The  first  is  modeled  by  a 
hydrodynamics  and  heat  flow  code  discussed  later  on.  This  code  outputs 
its  results  at  discrete  time  steps.  A  po s t-processor  code  models  the 
transport  of  radiation  through  the  plasma. 

The  time  scale  of  atomic  processes  is  very  short  compared  to  the 
time  step  of  the  hydrocode.  Ve  can  therefore  take  the  view  that  the 
atomic  processes  are  essentially  instantaneous,  or,  to  put  it  another 
way,  that  the  atomic  quantities  are  essentially  in  steady  state. 

In  this  view,  as  every  new  hydrocode  time  step  is  taken, 
ionization  and  excitation  processes  take  place  instantly  and  the  ionic 
species  densities  and  excited  level  populations  so  produced  remain  then 
unchanged  until  next  time  step. 

The  transport  of  radiation  through  a  plasma  and  the  ensuing 
emission  of  light  have  been  investigated  many  times  and  by  many  authors, 
and  several  different  approximations  have  been  made  at  different  times. 

Particularly,  a  point  that  does  not  seem  to  have  received 
overwhelming  attention  is  the  feedback  effect  that  the  radiation  field 
intensity  has  on  the  excited  state  populations  and,  thereby,  on  the 
various  parameters  characterizing  the  transport  (i.e.  emissivity  and 


3 


Absorptivity).  The  problem  here  is  that  the  population  of  the  excited 
levels  determines  the  light  intensity,  and  the  light  intensity  plays  a 
role  in  determining  the  population  of  the  excited  levels.  A  one-step, 
exact  solution  of  the  problem  is  therefore  extremely  difficult  in  the 
simplest  cases  and  impossible  generally. 

In  fact,  much  vork  has  been  done  assigning  some  guessed 
instantaneous  limiting  form  of  equilibrium  to  the  populations  then 
calculating  the  radiation  transport  and  emission  sustained  by  such 
populations.  The  criteria  proposed  for  the  guesswork  have  been  several, 
each  one  appropriate  to  a  particular  range  of  conditions.  Three  in 
particular  have  come  up  more  frequently:  LTE,  Coronal  equilibrium  and 
Escape  factors. 

The  first  two  imply  some  guess  as  to  what  the  radiation 
intensity  will  eventually  be,  approaching  black  body  in  the  first  case 
and  approaching  zero  in  the  second.  The  third,  the  Escape  factors 
method,  is  a  more  sophisticated  procedure,  first  developed  by  Holstein 
[1,2]  who  actually  solved  the  self-consistent  problem.  This  last 
method,  unlike  the  other  two,  takes  into  account  the  effect  of  the 
actual  radiation  field  on  the  excited  state  population.  In  it,  one 
first  calculates  the  coronal  population,  then  the  optical  depth  of  the 
plasma  in  the  center  of  the  line  with  such  population,  and  from  the 
depth  a  coefficient,  greater  than  one,  by  which  to  multiply  the  coronal 
value  to  obtain  the  population  of  the  excited  levels. 

But  in  fact  the  problem  can  be  exactly  solved  only  with  an 
iterative  procedure  in  which  an  initial  guess  on  either  the  population 
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or  the  light  intensity  is  followed  by  repeated  calculations  of  both  to 
follow  their  mutual  relaxation  until  some  convergence  criterion  is  met. 

To  accomplish  this,  it  is  necessary  to  know  the  emission  and 
absorption  profiles.  As  is  well  known,  atoms  do  not  emit  sharp,  delta- 
function-like  lines.  Several  phenomena  combine  to  distribute  the 
emission  (and  the  absorption)  over  a  profile  of  some  width.  The  three 
main  such  phenomena  are  natural.  Stark,  and  Doppler  broadening.  The 
natural  line  width  is  estimated  from 

AE*At~"h  where  At  ~  1/ASX 

AE~-fc  •  Aal~  10"*eV. 

As  for  the  Stark  broadening  [8,6]: 

AE~8.8  10“**  n9  eV 

nB  <  10ias“*  for  Hea  line  in  Al,  therefore 

AE  <  10“*  eV. 

The  Doppler  broadening  is  significantly  greater  than  either  of 
the  above.  Furthermore,  an  effect  that  has  not  received  much  attention 
in  previous  work  is  the  frequency  shift  due  to  the  expansion  velocity  of 
the  plasma.  Neglecting  it  can  be  acceptable  in  dealing  with 
astrophysical  sources,  but  is  clearly  inadequate  for  the  laboratory 
sources  of  x-ray,  in  which  we  are  primarily  interested,  since  sizeable 
plasma  flow  velocities,  in  units  of  local  sound  speed,  can  be  reached. 
It  is  a  general  property  of  ablative  flows  to  generate  velocities  from 
zero  to  supersonic,  and  the  difference  in  velocity  of  different  parts  of 
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the  plasma  is  such  that  the  difference  in  Doppler  shift  can  be  of  the 
same  order  of  the  local  Doppler  width.  This  entails  that  parts  of  the 
plasma  that  would  be  optically  thick  at  rest,  become  transparent  to 
radiation  from  other  parts  of  the  plasma  because  of  their  motion.  Ve 
have,  therefore,  used  a  multi-frequency  group  numerical  method  that 
makes  it  possible  to  include  Doppler  shift  from  flow  in  these 
calculations.  A  Doppler  profile  of  width  dependent  on  the  local  ionic 
temperature  is  assumed  in  the  present  work,  for  both  emission  and 
absorption  profiles,  assumed  identical. 

Since  the  code  transports  one  energy  group  at  a  time,  it  is 
convenient  to  take  a  somewhat  backward  point  of  view:  instead  of 
modifying  the  fequency  of  the  incoming  ray  to  the  frame  of  the  flowing 
plasma,  we  shift  the  absorption-emission  profile  of  the  moving  plasma, 
as  it  is  seen  in  the  lab  frame.  Since  the  flow  velocity  varies  with  the 
radius,  the  Doppler  shift  changes  through  any  given  shell.  It  is 
therefore  necessary  to  define  a  line  shape  that  is  an  average  throughout 
the  shell  (see  Appendix  A  for  the  derivation  of  this  average). 

Finally,  ve  have  dealt  with  the  angular  distribution  of  the 
intensity  with  a  novel  perspective:  in  spherical  symmetry  the  direction 
is  generally  described  in  local  spherical  polar  coordinates  in  which  the 
polar  axis  is,  at  every  point  in  space,  the  radius  through  that  point. 
Due  to  the  evident  azimuthal  symmetry,  only  the  colatitude  needs  to  be 
taken  into  account,  and  the  parameter  most  widely  used  is  the  cosine  of 
the  colatitude,  that  is,  the  cosine  of  the  sngle  that,  at  the  given 
point,  the  direction  considered  makes  with  the  radius;  this  cosine  is 
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generally*  end  will  be  throughout  this  work,  indicated  with  p  (see 
Figure  1). 

What  is  generally  done  in  numerical  transport  schemes  is 
discretize  p  into  'angular  groups.'  Instead  o£  introducing  angular 
groups,  we  follow  a  number  of  rays,  each  one  representative  of  a  group 
of  rays  and  characterized  by  its  distance  of  closest  approach  to  the 
center  of  the  sphere  (its  'impact  parameter')  in  its  path  through  the 
plasma.  Those  are  the  groups  of  characteristics  in  spherical  geometry. 
By  using  them,  we  eliminate  some  diffusion  truncation  error  that  would 
be  introduced  by  the  use  of  angle  groups. 

Summarizing,  the  novelty  of  our  approach  to  the  radiation 
transfer  resides  in  three  points  which  we  will  call  self-consistency, 
Doppler  effect  and  impact  parameter  groups.  We  shall  now  proceed  to 
discuss  the  numerical  algorithm  developed  to  model  this  physical 
situation. 

To  outline  the  procedure  briefly  (See  Figure  2): 

1)  The  atomic  properties  of  the  material  are  read  in.  The 
laser  and  target  parameters  are  read  by  the  hydrocode.  A 
hydrocode  run  is  made  separately. 

2)  The  hydrocode  output  is  read,  at  the  current  time  step,  in  a 
post-processing  mode.  Cold  zones  are  recognized. 

3)  Ionic  species  densities  and  electronic  density  in  every 
'hot'  zone  are  calculated. 

4)  The  transition  rates  are  computed.  The  population  is 
initialized  as  corona). 
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5)  Iteration  is  started:  radiation  flux  is  computed  and  then 
fed  into  the  rate  equation  to  obtain  the  next  iterate  for 
the  population*  then  the  radiation  flux  is  recalculated  and 
so  on  until  convergence  is  reached. 

6)  Back  to  point  2)  for  the  next  hydrocode  time  step. 


CHAPTER  2 


THE  HYDRODYNAMICS 

The  hydrodynamics  is  modeled  by  Chemise,  e  code  developed  by  R. 
Malone  to  model  the  implosion  of  targets  and  used  by  K.  Matzen  and  R. 
Morse  to  model  ablative  flovs.  The  code  follows  the  expansion  of  the 
plasma  in  time  dne  to  the  laser  energy  deposition. 

The  sphere  is  divided  into  n  shells,  or  zones,  aronnd  a  cold 
central  ball.  The  thermodynamical  quantities  are  assumed  constant 
throughout  any  given  shell. 

The  code  outputs  at  discrete  time  steps: 

Outer  radius  rj,  and  thickness  Ar^ 

Electronic  and  ionic  temperatures  Te^  and  Ti^ 

Ionic  density  Ni^  and  expansion  velocity  vx^ 
this  for  every  zone  k. 

This  output  is  read  by  the  transfer  code. 

A  zone  which  does  not  meet  given  requirements  on  electronic 
temperature  and  ionic  density  is  labeled  'too  cold'  and  considered 
inactive.  This  leaves  n'  'active'  zones.  For  every  active  zone  the 
densities  of  the  various  ionic  species  and  the  electronic  density  are 
calculated  using  the  ionization  and  recombination  coefficients 
calculated  by  Summers  [4,5].  In  the  treatment  of  the  excitation  states 
we  will  assume  that  ionization  and  recombination  play  a  negligible  role 
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in  swapping  electrons  fron  s  level  to  another,  and  we  will  neglect  it. 
In  this  way,  deternining  the  ionic  composition  of  the  plasma  and 
calculating  the  excitation  states  are  two  totally  decoupled  problems. 

In  the  first  round  of  calculations,  we  used  the  less  demanding 
(in  terms  of  computing  time  and  memory  requirements)  Two-level  atom' 
approximation.  This  well  known  simplifying  hypothesis  consists  in 
assuming  the  atom  capable  of  being  in  one  of  only  two  states:  the  ground 
level  and  the  (single)  excited  level  (henceforth  referred  to, 
respectively,  as  levels  1  and  2).  This  assumption  essentially  means 
transporting  one  line  at  a  time,  implying  that  the  radiation  field  due 
to  transitions  between  one  pair  of  levels  does  not  influence  the 
population  equilibrium  of  other  pairs.  Of  course,  this  assumption  is 
better  suited  to  situations  in  which  the  frequencies  of  the  transitions 
are  very  far  apart  and,  also,  the  ground  state  populations  is 
essentially  unchanged  by  the  excitation  state. 

Te  found,  as  expected,  that  it  is  not  quite  adequate  to  treat 
lines  that  are  close  enough  to  each  other.  Furthermore,  we  are 
interested  in  treating,  as  the  next  application  of  this  code,  x-ray 
lasing,  which  requires  considering  transitions  between  several  excited 
levels.  Ve  have  therefore  generalized  the  code  to  treat  multiple 
levels. 

In  the  following,  we  shall  refer  to  the  lower  and  upper  level  of 
any  one  transition  as  levels  1  and  2. 

The  information  on  electronic  and  ionic  densities  is  used  to 
calculate  collisional  excitation  and  deexcitation  coefficients  (ft**  and 


8*  i)  [6]  which,  together  with  Einstein's  coefficients  (A**,  Bis  snd 
Bji),  determine  the  rate  equations  for  the  atomic  transitions  between  an 
excited  level  and  the  groond  state. 
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where  $£  is  the  line  profile,  j£  is  the  direction-averaged  light 
intensity  and  n^  are  the  appropriate  densities. 

The  results  reported  here  are  fonnd  considering  only  transfer  of 
resonant  lines,  though  the  code  can  treat  transitions  between  excited 
levels.  Initial  runs  including  the  continuum  support  the  results  of 
experiments  being  conducted,  e.g.  by  H.  Epstein  [7],  and  show  that  the 
contribution  of  the  continuum  is  of  negligible  importance.*  The  excited 
population  is  obtained  by  assuming  steady  state  (i.e.  equating  the  right 
side  of  (II-l)  to  zero).  This  means  assuming  the  characteristic  time 
scale  of  the  hydrodynamic  processes  to  be  much  longer  than  the  time 
scale  of  atomic  transitions,  although  the  code  could  easily  treat  time 
dependence  of  level  populations  on  the  hydro  time  scale. 


*  See  Appendix  B 


CHAPTER  3 


THE  RADIATION  TRANSFER 

The  Bore  conventional  transport  schemes  working  in  spherical 
geometry,  generally  devised  to  describe  nentron  transport,  night  seen 
very  apt  to  model  the  transfer  of  light,  once  one  sees  the  light  as  a 
collection  or  beam  of  photon,  i.e.»  neutral  particles  only  interacting 
at  close  contact  with  scattering  centers,  and  in  this  respect  very  much 
like  nentrons.  What  such  schemes  generally  do  is  discretize  the  interval 
-l£|i£l  (where  p  is  the  cosine  of  the  angle  that  the  direction  considered 
makes  with  the  radius.)  into  a  number  of  angle  groups,  which  is  a 
natural  choice  in  planar  geometry.  But  in  curved  geometries,  a  straight 
trajectory  (such  as  an  unscattered  neutron  or  photon)  can  belong  to 
different  angle  groups  at  different  radii.  Therefore,  in  such 
geometries,  an  angular  group  treatment  requires  shifting  particles  from 
one  group  to  another  when  need  arises  due  to  the  change  in  spatial 
coordinate,  with  resulting  numerical  diffusion.  A  very  classical  example 
could  be  the  Sn  approximation.  But  taking  the  view  of  geometric  optics, 
which  is  how  we  generally  look  at  light,  it  seemed  more  appropriate  to 
do  the  contrary,  that  is,  adjust  the  groups  to  the  rsdial  position  in 
such  a  way  that  a  ray  (which  is  straight  in  the  absence  of  significant 
refraction)  remains  in  the  angular  group  throughout  the  transport. 

Having  chosen  this  approach,  it  seems  a  natural  choice  to 
characterize  the  groups  of  rays  by  their  'impact  parameter',  that  is. 
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the  distance  of  closest  approach  of  the  rays  representative  of  the  I 

groups  to  the  center  of  the  sphere.  This  procedure  has  the  fringe 
benefit  of  readily  yielding  the  intensity  as  a  function  of  radius  in  the 
projection  of  the  spherical  inage  on  a  plane.  That  is.  the  emitting  j 

sphere  as  seen  by  a  pin-hole  camera. 

As  for  the  choice  of  the  groups,  which  is  to  say  of  the  impact 
parameter  of  the  characteristics,  it  is  determined  largely  by  the  j 

hydrodynamics  description  used.  Because  our  plasma  is  modeled  as  a 
number  of  concentric  spherical  shells,  it  is  natural  to  assign  an  impact 
parameter  to  every  shell.  Our  choice  has  been  to  take  as  impact  j 

parameter  group  boundaries  the  inner  and  outer  radius  of  all  active 
shells,*  and  as  impact  parameter  of  the  representative  ray  the  average 
between  them.  This  leads  to  two  possible  situations,  as  illustrated  by  j 

Figures  (3  A  and  B):  every  ray  crosses  its  'critical'  rone,  that  is,  the 
zone  whose  inner  and  outer  radii  encompass  the  impact  parameter  of  the 
ray,  in  a  way  much  different  from  the  way  it  crosses  all  the  remaining  I 

zones  external  to  the  critical  one.  It  is  therefore  necessary  to 
specify  together  with  the  'normal'  transport  algorithm,  its  limiting 
form  for  the  special  case  of  the  critical  zone.  I 

As  for  the  boundary  conditions,  they  stem  from  the  requirement 
of  continuity  of  the  flux  at  every  boundary  between  zones. 


*  Extra  impact  parameters  through  the  cold  inner  plasma  can  be  added  at 
will  for  better  resolution  near  the  axis. 
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Particularly,  the  ingoing  flan  at  the  ooteraoat  surface  is  nnl. 
A  somewhat  less  obvious  condition  on  the  innermost  surface  is  that  the 
flux  from  inside  the  surface  toward  the  inside  of  the  first  shell  is 
nnl. 

So,  if  we  follow  a  'line  of  sight',  or  ray  (see  Figure  3C),  we 
start  with  zero  intensity  at  the  external  surface  and  build  up  intensity 
proceeding  inward  through  the  'normal'  zones,  then  through  the 
'critical*  zone  in  which  the  direction  changes  from  inward  to  outward, 
then  again  (but  this  time  proceeding  outward)  through  all  the  normal 
zones  and  finally  get  the  output  at  the  outermost  surface. 


Figure  3C:  Transport  Scheme 


CHAPTER  4 

THE  GEOMETRY 

Let  us  examine  the  geometry  in  detail.  Ve  have  n'*n-ncold 
zones,  and  we  set  n'+l  impact  parameters,  that  is,  through  the  n'  active 
zones  and  through  the  origin. 

For  every  group  l,  a  critical  zone  remains  defined  with  zone 

number  Rcrit^^  such  that  r(kcr ±t(l )-l)  <  dist  (Si)  <r(kcrit(j.))» 
dist (£ )  is  the  impact  parameter  of  the  I'tb.  line  of  sight.  For  purpose 
of  computation  we  will  consider  ngllg*2(n'+l)  angular  groups,  that  is, 
for  every  ray  we  consider  two  groups:  outward  (£)  and  inward  (narsj~*-+l)> 
so  that  a  group  turns  into  its  opposite  crossing  the  critical  zone. 

A  p is  therefore  defined  as  the  cosine  of  the  angle  between 
the  2,’th  ray  and  the  radius  at  the  surface  k.  An  angular  span 
is  associated  with  every  such  that 

UV-l)  Vk 

The  exact  prescription  may  vary,  as  long  as  (IV-1)  is  satisfied. 
Our  choice  has  been  as  follows. 

For  every  group  %,  let  us  define  two  more  directions,  one  with 
impact  parameter  r(kcr^^  (X)“l)  and  one  with  impact  parameter 
r(kCritU))»  that  i*  to  say,  ona  tangent  to  the  inner  surface  of  the 
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critical  zone  (for  the  £th  group)  and  one  tangent  to  its  outer  surface. 

Call  those  lines  ^  and  £  .  (See  Figure  4). 

<  > 

Evidently,  dist(^  )  <  dist(£)  <  dist(£  ).  Such  rays,  in  the  very 
sane  way  as  ray  l,  define  cosines  M^f  ^  and  ^  for  every  zone  k  and 
group  1. 

We  define  •* 

(IY-2) 

A»*k,£  "  V-i.l  “  ^k,A  Vk»£ 

It  is  apparent  that  m  and  Am  change  through  a  given  zone,  for 
any  given  ray.  It  shall  be  necessary  to  define  an  average  AMk^£,  to  be 
used  in  the  finite  difference  equation.  Again  the  main  requirement  is 
that  Ap^  i  obey  (IV-1).  For  the  sake  of  siaplicity  we  have  set  Am^  £  • 
Afk.jl*  tbat  is,  the  average  in  the  shell  equal  to  the  value  on  the  outer 
boundary.  This  choice  can  be  justified  as  follows:  for  any  ray  £  in 
the  critical  zone  the  external  surface  serves  both  as  entrance  and  exit 
surface.  In  all  other  zones  the  change  in  span  Am  through  the  zone  is 
limited. 

We  are  going  to  use  the  following  notations  for  intensities: 
out 

■  xE(xout,k  >0 
in 

^.k,^  "  xE(Rin,k  ’l) 

XE, k,  l  m  •▼•rage  of  lE(r,i.)  in  [r^-i.r^J 


Roy  (l)  is  representative  of  all  trajectories  with 

impact  parameter  comprised  between  drst  (i * )  and  dist(i>). 


Figure  4:  Impact  Parameter  Group  Boundaries 


CHAPTER  5 


RESULTS 

In  recent  jeer*  there  has  been  an  increasing  interest  in  intense 
beans  of  soft  x-rays  for  various  applications.  Very  bright  x-ray 
sources  are  needed,  for  instanoe,  for  the  pnrpose  of  etching 
semiconductors  now  that  a  reduction  in  size  of  one  or  tiro  orders  of 
magnitude  is  sought,  and  the  commonly  used  ultraviolet  beams  lack  the 
resolution  needed  for  such  miniaturization.  Laser  produced  plasmas  are 
a  prime  source  of  such  x-rays,  due  to  their  very  intense  brightness  and 
high  conversion  efficiency. 

One  apparatus  to  produce  such  x-ray  beams  consists  of  a  powerful 
Nd-glass  laser  whose  beam  is  directed  against  an  aluminum  flat 
target[7].  The  ensuing  short-lived  plasma  emits  very  strongly  in  the 
1-2  KeV  region  of  the  spectrum.  A  high  repetition  rate  makes  up  for  the 
short  life  of  the  plasma.  The  actual  target  consists  of  an  aluminum 
cylinder  that  rotates  and  advances,  not  unlike  in  an  old  phonograph,  and 
the  laser  beam  is  focused  on  a  40  |im  spot.  About  10  pulses  per  second 
are  shot  at  the  cylinder  surfaoe,  which  meanwhile  rotates  and 
translates  so  as  to  offer  a  fresh  target  to  every  pulse.  Each  laser 
shot  is  actually  composed  of  two  gaussian  pulses  (see  Figure  5),  a 
preheating  pulse  10  ns  long  and  depositing  1.5J  of  energy,  followed. 
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with  a  5  as  delay  between  peeks,  by  s  aain  pulse,  1  ns  long  end 
depositing  15J  of  energy.  The  long  preheating  pulse  creates  a  plasaa 
with  reduced  density  gradients,  thus  enhancing  the  absorption  of  the 
aain  pulse. 

The  outcoae  of  the  ezperiaent  is  that  aost  of  the  energy  in  the 
range  of  interest  is  eaitted  in  the  He-like  end  H-like  a-lines,  inasmuch 
as  the  continuua  is  negligible  and  the  higher  lines  are  at  least  one 
order  of  aagnitude  smaller;  the  line  ratio  of  Hea  to  Ha.  that  is.  the 
ratio  of  the  total  energies  eaitted  in  these  two  lines,  is  about  2,  and 
the  total  efficiency,  that  is.  the  ratio  of  total  eaitted  x-ray  energy 
to  total  deposited  laser  energy,  is  found  to  be  of  a  couple  of  percent. 

In  order  to  coapare  our  calculations  to  these  results,  we  aust 
first  aatch  our  input  parameter  set  to  the  experimental  conditions. 
Because  our  oode  gives  power  output  at  each  tine  step,  whereas  the 
experiaent  yields  total  energy  output  in  a  given  x-ray  line,  we  aust 
perform  at  times  integration.  A  siaple  procedure  to  do  this  is  to 
perform  a  rectangular  integration  over  all  the  time  steps  offered  by  the 
hydrocode.  This  being  extremely  expensive  it  seeaed  fit,  and  well 
within  the  precision  level  of  the  rest  of  the  approximations,  to  use 
only  a  limited  number  of  representative  steps. 

The  hydrocode  being  adaptive,  the  time  step  is  automatically 
adjusted  so  that  an  approximately  constant  change  in  thermodynamical 
properies  occurs  from  step  to  step.  It  seems  therefore  reasonable  to 
choose  a  constant  pace  in  selecting  representative  steps,  like  every 
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100th  or  overy  300th  *tep.  The  actual  number  can  be  selected  by  the 
user  to  fit  best  the  problem  they  are  treating. 

To  match  the  code  to  the  experimental  conditions  requires  some 
more  careful  justif i cation.  In  the  first  place,  our  code  is  designed  to 
model  spherical  targets,  whereas  the  experiment  described  above  involves 
flat  targets.  It  has  been  shown,  however,  [10]  that  the  hydrodynamics 
of  the  expansion  of  a  laser-heated  flat  target,  having  a  spot  diameter 
D,  can  be  modeled  with  good  accuracy  as  a  conical  section  of  a  spherical 
expanding  plasma  (see  Figure  6). 

More  precisely,  starting  with  a  spherical  target  with  radius  R«D 
and  heating  it  with  a  spherically  symmetric  laser,  with  the  same 
intensity  as  for  the  flat  target,  the  portion  of  the  ensuing  plasma 
comprised  inside  a  cone  of  semiaperture  30°  and  vertex  in  the  center  of 
the  sphere  will  model  very  closely  the  actual  plasma  generated  by  the 
flat  target.  This  model  is  particularly  good  near  the  axis  of  the  cone, 
whereas  it  fails  at  the  edge  of  it,  as  one  would  reasonably  expect.  Of 
course,  the  laser  energy  deposited  on  the  whole  model  sphere  is  sixteen 
times  the  energy  actually  deposited  on  the  real  flat  target,  and  this 
must  be  taken  into  account  in  choosing  input  power  parameters  for  the 
code. 

The  absorption  routine  in  the  hydrocode  considers  only  inverse 
bremstrahlung,  and  actually  deposits  only  some  fraction  of  the  total 
laser  energy.  To  fit  Epstein's  experiments  best,  we  adjusted  the  total 
laser  power  so  that  the  predicted  energy  absorption  would  equal  the 
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absorption  actually  measured  in  the  experiment.  In  Figure  7  we  can  see 
what  the  energy  deposition  looks  like  as  a  function  of  radios,  around 
the  peak  of  the  laser  pulse.  The  cut-off  at  the  critical  surface  is 
▼ery  evident.  Thanks  to  preheating,  a  large  fraction  of  the  laser 
energy  (58%)  gets  deposited.  Under  such  conditions  neglect  of 
absorptive  nechanisas  other  than  IB  is  believed  to  be  a  reasonable 
approximation. 

Table  1  shows  the  paraaeters  of  the  run.  Figures  8  to  11  show 
the  spatial  behaviour  of  the  plasaa,  as  calculated  by  the  hydrodynamics 
and  heat  flow  code.  Figure  8  shows  the  electronic  temperature:  it 
follows  somewhat  the  shape  of  the  laser  deposition  around  the  critical 
surface.  Then,  proceding  outward,  one  may  see  that  the  high  electronic 
thermal  conductivity  tends  to  homogenize  the  temperature. 

In  Figure  9  we  show  the  ionic  temperature.  Ions  are  heated  by 
cross-relaxation  with  the  electrons,  and  by  viscosity  when  shock  waves 
pass  through.  As  an  example  of  viscous  heating,  there  is,  in  Figure  9, 
a  second  peak  at  about  0.0 6  cm  radius  that  clearly  cannot  be  the  result 
of  electron  heating.  This  suggests  that  a  shock  is  building  up.  This 
impression  is  confirmed  by  a  close  scrutiny  of  Figure  10.  A  small 
density  peak  is  building  up  at  the  same  radial  position,  because  of  the 
sudden  heating  by  the  short  main  pulse. 
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Target: 

Material 

Size 

Initial  temperature 
Laser  Pulse: 

Shape 

Energy 

Full  width  at 
Half  maxi  mini 

Delay  between  peaks 

Wavelength 

Total  absorbed 
energy 


TABLE  1 

Aluminum 
Radius  *  40  pm 
0.25  eV 

2  gaussians 

El  -  2.6  J 
E2  -  26.0  J 

FWHM1  «  10.0  nsec 
FWHM2  «  1.0  nsec 

"  5.0  nsec 

*  1.06  pa 

-  16.5  J 
(-58%) 
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Particularly  noteworthy  is  the  expansion  velocity  (Figure  11) 
which  is  above  sonic  in  all  that  part  of  the  plasma  that  we  may  expect 
to  contribute  most  to  the  radiation,  that  is,  the  part  that  is  dense  and 
hot  at  the  same  time.  This  points  to  the  need  for  a  treatment  of  the 
Doppler  shift.  Ve  can  see  that  this  radiating  part  is  a  relatively  thin 
shell  around  the  critical  surfaoe,  since  density  drops  sharply  with 
radius.  The  effeot  of  the  two  pulses  is  very  clear.  The  velocity  at 
large  radius  is  smaller  than  the  velocity  close  to  the  oritioal  surface: 
we  must  bear  in  mind  that  the  three  curves  shown  are  at  three  times  very 
close  to  each  other,  around  the  peak  of  the  main  pulse.  Therefore  the 
effects  of  this  expansion  have  not  reached  the  more  peripheral  plasma, 
whose  expansion  velocity  is  still  the  one  due  to  the  preheating  pulse. 

This  graph  also  supplies  evidence  for  the  explanation  of  the 
second  ionio  temperature  peak:  following  the  velocity  profile  inward, 
the  velocity  suddenly  steepens  around  r»0.06  cm.  A  faster  expanding 
plasma  is  pushing  against  a  slower  one,  and  a  shock  wave  is  to  be 
expected. 

Figure  12  axe  the  first  results  of  the  post-processor.  As 
mentioned  earlier,  the  code  recognizes  'cold'  zones,  and  totally 
neglects  them.  Therefore  no  electronic  or  ionic  densities  are 
calculated  fox  such  zones.  That  is  why,  in  Figures  12,  such  densities 
start  abruptly  with  a  very  high  value  near  the  critical  surface.  Ve 
will  return  to  these  graphs  later. 


12B 


41 


In  the  present  calculation  wo  haw*  aaed  100  photon  energy  groups 
per  line*  200  for  the  doublet.  Although  to  prediot  the  total  energy 
output  fewer  groups  would  suffice*  we  found  that  a  much  smaller  number 
of  them  would  give  a  wary  rough  definition  to  the  line  profiles. 

.the  code  checks  itself  for  consistency,  comparing  the  net 
collisional  pumping  power  in  every  zone  with  the  net  flux  of  radiation 
out  of  the  zone. 

As  it  turns  out*  the  radiation  output  is  perfectly  consistent 
with  the  collisional  power,  zone  by  zone  and  globally.  It  is*  in  fact* 
a  necessary  condition,  given  this  physical  situation,  that  net  radiation 
output  and  net  pumping  power,  for  every  zone*  be  equal  to  each  other. 

On  the  other  hand,  if  we  look  at  the  laser  energy  deposition 
(Figure  7).  we  notice  that  almost  all  the  energy  is  deposited  in  a 
region  comprising  only  a  couple  of  zones.  This  seems  to  be  at  the  very 
lowest  boundary  of  acceptable  resolution,  even  taking  into  account  the 
fact  that  the  various  zones,  passing  through  this  'deposition  region'  as 
time  passes,  'average'  in  some  sense  the  state  of  the  region.  It  is 
necessary  to  have  a  finer  grid,  so  that  the  deposition  is  smeared 
through  a  large  number  of  zones,  and  w&  plan  on  doing  so  in  the  very 
near  future. 

Figure  15  shows  the  emitted  spectrum,  integrated  over  time.  An 
integration  over  photon  energy  shows  that  Hea  and  Ha  lines  account  for 
more  than  90%  of  the  total  emitted  energy,  the  line  ratio  of  Hea  to  Ha 
is  2.5-3.00  and  the  total  efficiency  is  4.5-5%:  the  uncertainty  in 
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those  figures  steals  fro*  the  fsot  thst  those  quantities  are  rather 
sensitive  to  the  exact  aaionnt  of  energy  deposited  by  the  laser,  and  the 
absolute  value  of  this  deposited  energy  is  not  easy  to  measure  in 
experiments.  Those  results  are  in  good  agreement  with  Epstein's 
experiments,  especially  the  line  ratio,  which  is  a  very  sensitive 
measure.  Details  of  the  time  integrated  line  shapes  are  given  in  Figure 
Id  to  18. 

Figures  19  to  24  are  in  the  same  vein,  but  referring  to  the 
instantaneous  emission  at  peak  power.  The  narrow  shape  of  He-like 
lines  (compared  to  H-like,  that  is)  is  explained  with  the  help  of  Figure 
12B.  As  we  see,  He-like  ion  density  is  very  peaked,  and  rolls  off  very 
quickly.  H-like  ion  density,  on  the  contrary,  increases  at  first  with 
radius,  and  then  rolls  off  more  gently.  Te  can  hence  expect  that  the 
He-like  emitting  zone  would  be  deeper  inside  the  plasma  and  extend  over 
a  smaller  thickness  of  plasma  than  the  H-like.  The  average  flow 
velocity  of  the  He— like  emitting  plasma  is  therefore  mar-  “dly  smaller 
than  the  average  flow  velocity  of  the  H-like  emitting  one,  this 
naturally  implies  a  smaller  Doppler  shift  and  hence  less  smearing  of  the 
line  shape. 

To  understand  the  result  of  the  multi-level  calculation  of  the 
Ha  doublet,  we  have  to  look  at  Figures  25.  Those  line  shapes  have  been 
obtained  by  arbitrarily  reducing  the  expansion  velocity  of  the  plasma 
(to  0,  25%,  50%  and  75%  respectively),  and  therefore  the  amount  of 
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Figure  16:  Time-integrated  He  Line 
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Doppler  shift.  As  shift  increases  the  second  line  of  the  doublet  gives 
up  its  energy  to  the  first  line. 

This  can  be  explained  as  follows:  light  eaitted  in  the  second 
line  at  soae  point  in  the  plasaa  is  reabsorbed  in  the  first  line  at 
points  further  out.  therefore  contributing  to  enhance  the  first  line's 
excited  level  population,  whereas  light  coaing  froa  inside  cannot  excite 
the  second  line  at  outer  points.  The  global  effect  of  this  phenoaenon 
is  to  shift  radiation  froa  the  second  to  the  first  line  in  the  doublet. 
This  is  aore  markedly  so  the  larger  the  Doppler  shift,  until,  if  the 
plasaa  is  expanding  fast  enough,  the  second  line  coapletely  disappears. 

Of  course,  the  opposite  is  true  for  light  transported  inward,  so 
that  we  can  think  of  the  population  of  the  excited  levels  of  the  two 
lines,  when  the  two  lines  are  treated  together,  as  'tilted'  toward  the 
second  line  in  the  inner  plasaa  and  toward  the  first  line  in  the  outer 
plasaa,  as  compared  to  the  populations  we  would  find  treating  one  line 
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CHAPTER  6 


CONCLUSIONS 

A  code  has  bean  developed  that  is  both  economical  and  accurate. 
It  is  designed  so  that  more  accuracy  can  be  obtained,  of  course  with  a 
higher  cost,  or  it  can  be  made  more  economical,  whenever  less  accuracy 
is  needed.  The  results  are  in  good  agreement  with  the  experimental  work 
to  the  extent  of  the  experimental  resolution.  Unfortunately  experiments 
are  as  yet  unable  to  resolve  line  shapes.  The  detail  with  which  the  line 
shapes  can  be  predicted  by  our  code  have  already  prompted 
experimentalists  to  try  and  improve  the  resolution  of  their  experiments 
to  measure  line  shapes. 

Whenever  experimental  data  on  line  shapes  will  be  available,  we 
will  check  our  results  against  it.  We  plan  also  on  using  our  model  to 
predict  diagnostics. 
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APPENDIX  A 


VARIOUS  DERIVATIONS 


a)  The  Rate  Equation 


From  (II-l),  with  the  condition 


dna 


dt 


(steady  state  assumption),  we  obtain 


/• 


JgdE)ni  *  (Ax 


»+Qiine+Bia  I  <| 

J  r\ 


^  gdE)n  2 


with  the  relations  (9] 


Si  Eh 

On  *  Du  -  exp  { - 

8*  KT« 


8s 

Ba,  * - B»a 

8a 


and  gj  are  the  statistical 
weights  of  levels  1  and  2) 


B,3 


-i  *  , 

h  c  Aa  a 


2  (Eaa)* 


Therefore  we  can  write  a  system  of  2  equations  in  nx,  m: 


Oisna  +• 


8s  f 

- Baa  / 

8a  J 


$  gdE 


Ea: 


Asa  +  Qisne - e 


8  a 
8s 


ETe 


/ 


+  B»i  I  ^  gdE 


i 


j 


a 
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m 
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ni  +  o»  * 

which  yields  nj.  and  n«. 
Also 


£  JEdE  ■ 


<fcm  J“  AEa 


wksre  the  stim  is  over  photon  energy  groups 


b)  The  Finite  Difference  Equation 

We  will  use  the  following  notations: 

Rlk  *  rk-i»  “  rk 

Alk,  A2^  and  Vj.,  area  of  the  inner  and  outer  surfaces  and 
voluae  of  the  shell. 


I 


J 


The  subscript  'in'  refers  to  the  surface  through  which  the  ray  enters 
the  zone,  'out*  to  the  other  surface. 
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The  energy  flax  Fg(r)  is  defined  as 

/ 


■/: 


a  Ig(r,Q)dQ 
AD 

Simple  energy-conservation  considerations  yield 


V->-  Source 

•  f£  *  -  - 

Unit  Volume 

Aaxna  E1X  ^g”  (Bjjn^B^nj)  EX1  ^  g  Jg 

where 

^xa,  Bxa,  Baa  ■  Einstein  coefficients 
<f>  g  *  Absorbtion  profile  ■  Emission  profile 


f  IE<r* 
J  AQ 


Q)  dQ 


Jw(i) 


f  dH 

J  AH 


Let  as  consider  a  spherical  shell,  with  constant  thermodynamical 


properties,  bounded  by  two  spherical  surfaces  of  radii  R1  and  R2. 


The  aet  flax  of  Fg(r)  oat  of  the  shell  is 


FE(r)  ‘  a  dA 


I  F(B2) 'r  dA  -  I  F(B1) *r  dA 
*^S2  SI 


2n  A2  f  |ia  Ig(B2,(ia)  d|ia  -  2n  A1  I  pia  Ig(Rl, |ii)d|ta 


/■ 


where  p  ■  n(r),  |ia  ■  |i(Rl),  |ia  *  n(K2). 


Oa  the  other  haad 


2jt ^d|i 


V  *  Fp  ■  Aaa  aa  Eaa  <p 


in 


~  (Baajia  -  Bal|ia)Eaa  0  g  - 


/ 


2n  /  Ip(r,|i)dn 


4jt 


Now: 
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2)i  A 2 fij  IgCHlf  )ij)  dfij  2n  A1  /  Hi  IE(Rl,|ii)  dpx  » 
[ff/  ^*in»E*i  (Bxanx-Ba  xna)  \ 

J \J  \~  *E  *  E,‘  iE  )  d“ 


2nj'jj'(s  tg  -  a  *E  IE(r, n>)  dn  | 


dV 


where 


Aa xnaEax 

S - , 

4rr 

BXanX“BaXna 

H  - - Eax 

4jt 

This  yields  the  finite  differences  equation. 

A2  iitfi  Ia  -  A1  px£  A|ix>£  Ix ji 

*  V  <P  E  (  SAjlj  -  H  A^£  l£  ) 

where  S  and  H  are  constants  and  A(i£,  I  are  properly  defined 
through  the  volume. 

Finally: 


average  s 


A2  Vijl  dpa^  Ia>£-Al  pXj£  Apx>£lx>£  =  V  «j>“  (&-H  1 1  )  A|iJl 


Tiro  cM.se a  are  possible: 
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1)  Forward,  that  is  |ix  and  |is  >  0 

Then  I is  specified  as  boundary  condition 

2)  Backward,  that  is  and  n*  <  0 

here  the  boundary  condition  is  on  I*  £ 
Accordingly,  we  have  the  following  two  equations: 

Al|ix  oAjii  o  V<J>®Ajtj i 

1)  I»  £  - -  Ilfs,+ - 

A2|ia>JiAjia  tl  A2|i»>£Ani>£ 


(s-aij) 


m  A2pa#JlA|ia  o 
2)  lt'2,  +"  ' 

Alpi  ^Ap!  ^  '  Aljix^Api t% 

(a*)  In  the  Critical  Zone 

For  the  limiting  case  of  the  critical  zone,  notice  that 

•H'  *  "  H  •  •  A|i£  . 

the  fozmula  used  is 


(Sr-BlZ) 


m  m 

I*,f  *  I»,l '  + 


V 


2S  -  H(l“  +  1 *  ) 

l  l 


b  * ')  Is  A  Cold  Zone 


m 


If  S  *  H  -  0,  then  I„  (r)  -  I 


R 


in 


in,£ 


therefore: 


m  ra 

l  *  ^n,  i 


AiaAR 


_m  ®in 
*out,£  “  *in,£  “T" 


R 


out 


J 

< 

I 


--I 

■ 


M 


M 


M 


M 


fl 


J 


jp 
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So: 


♦ 


£ 


[A  +  lUr-TjL)] 


[(A+B(r-ri)] 

dr  » 


1 


2AEQAr 


1 


2ABoAr 


1 


B 


1 

B 


2 


jr 


2 

& 


[A+B(r-rx)] 


A+BAr 


a 

d[A+B(r-rx)] 


2AEoBAr 


APPENDIX  B 


CONTRIBUTION  OF  THE  CONTINUUM 


*)  go&iiassa  gglajlaa  msA  absorption 

1)  Absorption 


KEcont  -  sEff  +  iE 


bf 


IEff  -  ck 


f  f 


W  — 

lE  J  ^N. 


Ckff  -  2.4234  10 


-37 


KEbf  -  Ck*>f 


N 

E3~  i 


(i~l)4 


Nj 


y 


Ck*>f  -  1.9924  10“14 


If  we  assume  ni<Q  s  »  Nj,.  n^j  -  0  y  j  ^  G.S.  then  the  above 
equation  simplifies  to 


Kgbf  -  Ckb* 


(i-l)4 

^i,  g.s.)5 
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Since  we  are  considering  an  interval  aronnd  ve  make  the 

further  simplification  that 


1  1 
13  Ei.j3 


in  both  formulas  for  and  Ke^* 


2)  Emission 


eEcont  -«Eff  +«Efb 


C£^  -  5.1965  10“14 
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Again,  assuming  ni#g>s<  “  Nj  ,  ni,j  ”  ®  Vj  #  6.S.,  and  neglecting 

ADj 

the  correction  for  the  lowering  of  the  continuum  (exp  [ - ]=1): 

Te 


The  results  show  an  energy  emitted  in  the  continuum  per  unit  eV 
that  is  much  smaller  than  the  one  in  the  line. 

Extrapolating  the  result  to  the  entire  range  1500  -  2500  eV 
gives  a  total  energy  emitted  in  the  continuum  in  such  interval  **  about 
100  times  smaller  than  the  energy  emitted  in  the  lines. 

Econt  "  9.5  10"7  =>~  9.5  10"4  J 

from  hv  -  1500  to  hv  «  2500. 

Eune  *  9.5  10“2  from  hv  -  1500  to  hv  -  2500. 

**  For  the  examples  we  treated. 
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